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Решение интегральных уравнений 
в системах компьютерной алгебры 


Представлены инструментальные средства для системы алгебраического программирования АР5, 
Маре, Мафетайса и других систем компьютерной алгебры. Этими средствами конструируют программы 
для таких систем. Эти программы имеют следующие вход и выход: 


ПМРОТ = (Е[›]= 0, [4,6], и), где Ву] = х, ув, |4’ Кбьв, У(0) 46, 
ОПТРОТ = у, = с с,х+ ...+ вх" = у(х ) = зоме( Е у]=0), хе [а,6] , 
|| УС®) — у»(®) [авы / Ме о, с в | У(Х) - (0+... + сх) | саы = ОЦ). 


Этими инструментальными средствами сконструирована процедура системы АР$. Эта процедура 
решает линейные интегральные уравнения вида 


Ау+Цу] +8=0, где А = зоуе(Ву + Ну] +/=0), Ду] = Киу] +...+ КУ], 
НЫ] = К» + ... + КОЛ, КОЛ = 59“ Кд д 4, 
ядра К{х,)) и коэффициенты сдх), 4х), © , В, Г- многочлены. Доказана эффективность этой процедуры. 


Доказана эффективность метода — основания этой процедуры. Процедура доказывает эффективность 
этих средств. 


Введение 


Задача. Построить инструментальные средства для системы алгебраического 
программирования АР$, Маре, Мафетайса и других систем символьных вычис- 
лений на компьютерах — систем компьютерной алгебры (СКА). Этими средствами 
конструируют программы для таких систем. Эти программы имеют следующие вход 
и выход: 


ПУРОТ = (51 = 0, [4,5], я), где Е] = Ах, 69, [КС в 6) 40, (1) 
ОПТРОТ = у, = с + сх +... сих" = у(х) = з9уе(Е[ у] = 0),х Е [а,ё]. (2) 


Коэффициент оптимальности алгоритма, реализованного в такой программе, на 
интегральном уравнении (1) в пространстве Х = Са] ограничен 


С‚(а|еогит, Ту] = 0, Х) = [у - ух, Е о, де п|| У) = (С +...+ сих”) |5 = ОС). (3) 


Актуальность задачи. Маре, Мафетайса и другие СКА вычисляют решение 
обыкновенных дифференциальных уравнений в аналитическом виде. 

СКА не решают аналитически функциональные уравнения других типов. 

Интегральные уравнения моделируют явления окружающего нас мира более 
широкого класса, чем обыкновенные дифференциальные уравнения. 

Возможности классического математического обеспечения. Программно 
реализованы, как правило, методы с насыщением. По этим программам вычисляют 
сеточную аппроксимацию {у; = у(х/} 0” решения у(х) интегрального уравнения. Для 
исследования этого сеточного решения {у; = у(х;)} о’ на отрезке [а,Б] его интерпо- 
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лируют — преобразуют в алгебраический многочлен Р„[{уй =о']. Если функция у — не 
целая и узлы сетки {хи о” — равноотстоящие на отрезке [а,5], то с ростом и алгебраи- 
ческий многочлен Р;[{у(х)} 0] (интерполирующий функцию у в этих узлах) часто не 
сходится к функции у в пространстве Сраь. 

Методы ряда Чебышева решают интегральные уравнения [1, с. 275-348] и вы- 
числяют аппроксимацию частной суммы ряда Фурье — Чебышева решения. СКА не вы- 
полняют преобразования этих методов. 


1 Инструментальные средства 


1. Операторы на языке АРГАМ системы алгебраического программирования АР5З. 
Эти операторы выполняют преобразования а-метода Дзядыка [2]. 
2. Процедуры из этих операторов для решения линейных интегральных уравнений 


Ау Пу] +8=0, (4) 
где интегральный оператор имеет вид 
ЦУ=КиУ +... + К, КУ = [ку © ККх,0 У а (5) 


ядра КХх,6), коэффициенты сх), 4(х), А и свободный член © — многочлены (линейные 
интегральные уравнения с многочленными коэффициентами — ЛИУМК) следующих 
типов: 

— без особенности на отрезке аппроксимации [а,в] [3], 

— с одной регулярной особой точкой ноль. 

3. Метод конструирования процедур для решения интегральных уравнений: 

— аппроксимация исходного уравнения последовательностью ЛИУМК (4), 

— решение ЛИУМК этой последовательности по указанной выше процедуре. 


2 Решение интегральных уравнений типа А5 
ПМРОТ = (Ау + Цу] + & =0, А = зауе(В у + Ну] + }= 0), [а,6], п), (6) 


где ©, /, В — многочлены, Гу], Н[у| — операторы вида (5). 
Метод 1. Распространение а-метода Дзядыка решения ЛИУМК Вольтерра [2]. 
1. По методу [3] преобразовать ЛИУМК В у + НуУ]| + Г = 0 системы (6) в 
многочлен 


А„ = те®о4 [31(В у + Ну] + Т= 0, [а,Ь] п). (7) 
2. Преобразовать уравнение А у + Ду] +2 = 0 системы (6) вЛИУМК 
Ау + Пу] +8=0. (8) 
3. По методу [3] преобразовать ЛИУМК (8) в многочлен 
у = шефо4 (А, у + Ду] +8 =0, [а,Б] п). (9) 


Алгоритм 1. 
1. По алгоритму [3] преобразовать ЛИУМК Ву + Н| у] += 0 (6) в многочлен 


А„ = а ог из (В у + Шу] + 1= 0, [а,6], п). (10) 
2. Преобразовать уравнение А у + [у] +2 = 0 системы (6) вЛИУМК 
Ау + Цу| +2=0. (11) 
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3. По алгоритму [3] преобразовать ЛИУМК (11) в многочлен 
у„ = ао (Ан у + Ду] +5 =0, [а,Ь], п). (12) 


Лемма 1. Метод 1 эквивалентен алгоритму 1 — тождественны многочлены у» 
(9) и (12). 

Доказательство. Согласно определению метода [3] и алгоритма [3], многочлен 
А, (10) тождественен многочлену А, (7). Поэтому тождественны уравнения (8) и (11). 

Следовательно, согласно определению метода [3] и алгоритма [3], многочлен у, (12) 
тождественен многочлену у, (9). 

Из спецификации алгоритма 1 и алгоритма [3] следует такое утверждение. 

Лемма 2. Все преобразования алгоритма 1 — алгебраические. 

Структура данных на входе АРГАМ-процедуры. 

1. ЛИУМК Ву + Ну| + {= 0 системы (6) имеет вид требуемый процедурой [3] 


ОМК := В *у + Е ор(с_ 1, а Т,К Т*у, + 
...+ ИИ ор(с 5,45, К 5*у, И +Г= 0, (13) 


где у, { — атомы. Эта процедура решает ЛИУМК. Ядра К 1....,К_$ — термы, имеют 
естественный вид и аргументы х ‚ {— атомы. 

Коэффициент В, сводный член © и пределы интегрирования с_1....,с_$, 4 1,..., 
4 $ - термы, имеют естественный вид и аргумент х - атом. 

2. Уравнение А у + Ду] + & = 0 системы (6) определяет оператор 


Гу: = тЕ ор(с_1, 41, К 1 *у, 1 +...+ тЕ ор(с р, ар, Кр *у, И +в. 


3. Отрезок аппроксимации [а, 6] определяет список (а, Б). 
4. Параметр и является целым числом. 
АРГАМ-процедура 1. Алгебраическая спецификация алгоритма 1. 


А п := уе 1 а(ГЛОМК ицегуарп); /* А_п */ 
ГОМК 1 := (А п*у-+Гу=0); /* Ап*у+ [у] + #=0*/ 
у п := зауе 1 ч(ГТОМК _1Тимегуа|,п); /* уп */ 


Замечание 1. АРГАМ-процедура зо]уе_1 и — реализация алгоритма [3]. 
Структура выхода операторов АРГАМ-процедуры. 
1. Многочлен А_п (10) имеет коэффициенты - числа и вид 


о РЕВ, (14) 


2. Уравнение ГТОМК 1 (11) имеет вид (13). 
3. Многочлен у_п (12) имеет коэффициенты -— числа и вид (14). 


3 Эффективность АРГАМ-процедуры 
Теорема 1. Для сложности АРГАМ-процедуры справедливо тождество 


О(аеогийт_1(А у + Ду] + 2 = 0, А = зоуе(В у + Шу] +Х= 0), [а,Б], п)) = 
О(Т) + О(зауе_ т и (Ву+ Ну] +Х=0, [а,6], п)) + 
О(зо1уе_1 и (Алу Му] +5 =0, [а,6], п)) , 


где О( зауе 1 и (Ву-+ Ну| + Г = 0, [а,6], п)) — (полиномиальная) сложность решения 
ЛИУМК Ву + Ну] += 0 процедурой за уе_1 и [3]. 

Доказательство. АРГАМ-процедура 1 — линейная. Сложность линейной про- 
цедуры тождественна сумме сложности её частей. Сложность формирования ЛИУМК (8) 
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не зависит от параметра и. Согласно оценке [3], процедура зойуе_1 и имеет полиноми- 
альную сложность. 

Теорема 2. Пусть числа — константы операндов входа АРГАМ-процедуры 
являются целыми или рациональными. 

Тогда АРГАМ№-процедура выполняет преобразования алгоритма [ точно. 

Доказательство. Согласно лемме 2, преобразования алгоритма 1 — только алгеб- 
раические. При алгебраических преобразованиях система АР$ (и другие СКА) выпол- 
няет арифметические операции с целыми и рациональными числами (константами 
операндов этих преобразований) точно. 

Модельное интегральное уравнение метода ряда Чебышева [1 с. 305] 


Ау- |" (-Вуфанх=0, (15) 


А = зоуе(у + 4 [© (х-0 у @- 112-12 =0)= с0з(х)?/2. (16) 
Описание системы уравнений (15), (16) на языке АРГАМ. 
ргосезз[1] := (ГТОМК := (у + (-1/2) хх^2 + (-1/2) + 
шШЕ ор (0, х, (4х (х+ -1х 9) ху, © =0); 
Ку := шЕ ор(0, х, (-1 х (х+ -1х1)) хху, © + -1 хх; ...); 
Результаты решения системы (15), (16) процедурой. 
п := 3; имегуа[ := (-1, 1); 
А_п = гаК-6,17) хх ^ 2 +гак3З3,68) ; 
ГОМК 1 = (гак-6,17) хх^2 + гаКЗ3,68)) ху + 
шШЕ ор (0, х, (-1 х (х+-1х 9) ху, + -1 хх =0); 
у п = га{(57664,6271) хх ^ 3 + гак-2040,6271) хх; 


4 Эффективность алгоритма 1 


Построенная выше АРГАМ-процедура и теоремы 1, 2 доказывают следующие 
утверждения. 

Теорема 3. Сложность алгоритма 1[ — полиномиальная по параметру п. 

Теорема 4. Пусть числа — константы операндов входа алгоритма 1 являются 
целыми или рациональными. 

Тогда СКА выполняют преобразования этого алгоритма точно. 

Коэффициенты Фурье — Чебышева решения уравнения (15) 


У(о= зоуе( соз(х)?/2 у — [6 (х- Ву @и-х=0)=2 0 /\соз(х)’ (17) 
на отрезке [-1,1] принимают следующие значения а2 {(у, [-1,1]) = 0, 


{азна(» [По = $7,3, 2,6, 0,65, 0,14, 0,027, 0,0049, 
0,00084, 0,00014, 2,3 10°, 3,6 10°, 5,6 1077, 8,4 10°} (18) 


и регулярно убывают с ростом параметра 21+ 1 


Г |+ 3, [+1 [а +10», [ПП] = 4 
‘а но = {0,36, 0,25, 0,22, 0,19, 0,18, 0,17, 0,167, 0,164, 0,156, 0,155, 0,15}. 


Следовательно, справедливо тождество 


Е 0... сп || УСО - (С +...+ сих”) [сньи= [@2 ен 1, [-1,1]) [1 + В„), 
где В» < 1/(1 — Че-102), {ИА — ад} = {0,6,\ 0,3, ..., 0,17}. (19) 
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Норма в пространстве Х = С. погрешности аппроксимации многочленом у, = 
= аи (А у — [6 (х- В Уд &-х=0,у+4 [0 (<-0 У @- 2х — 12 = 0, [-1,1], п) 
решения у (17) уравнения (15) принимает следующие значения 

{Пу уэна |х= | У-— зна |х н-о = 46, 1,8, 0,25, 0,062, 0,0093, 0,0017}. (20) 

Из тождеств (18), (19) и (20) следуют заключения. 

Вывод 1. Главной частью коэффициента оптимальности (3) алгоритма [ на 
интегральном уравнении (15) в пространстве Х = С 1] является отношение 

[р уж Пень / | 42 чпу 1, [-1, 1) [= С + азще-ань 

где {С+ оон 0 = {2,6, 2,8, 1,8, 2,35 1,9, 2}. 

Вывод 2. Для коэффициента оптимальности (3) алгоритма 1 на уравнении (15) в 
пространстве С! - справедливо тождество 

С; аогит_1, едаайоп (15), С|-1,1]) = (С + ор-рни / (1+ В») , 


где В, < (1 — 9-12), {ИА - авео = {0,6 0.3, ..., 0,17}. 

Замечание 2. Минимальное значение коэффициента оптимальности алгоритма 
преобразования системы уравнений (6) в алгебраический многочлен - 1. 

Замечание 3. Коэффициент оптимальности метода ряда Чебышева на задаче 
Коши, эквивалентной системе (15), (16), принимает близкие значения [1, с. 306]. 


5 Вариант применения метода Галеркина 
ПМРОТ = (Ау + Цу| +8=0, А = зоуе(В у + Ну] + Х= 0), [а,6], п) , 
ООТРОТ = баегкт(А у + Цу]| + ® = 0, А = зоуе(В у + Н [| +Х= 0), [а,Б], п) = 
у = 0\е(5„[А„ (у Е Ни[а,Ь]) + Ы[ 0» е Нё“а,6]] + 8] = 0), 
А, = зоуе($[В (у, Е На,Ь]) + Н[@», е Нца, += 0). (21) 
где 5„[у| = а0(у, [а,6]) свеб(0,2(х)) -...+ а,(у», [а,6]) свеБ(и,2(х)) : [2(а,Б; р) >> На, 6] 


свеб( р, х ) = соз(1 агссо$( х)), 2(х) = 2 (х-а) / (Ь-а) — 1, 
(и Е Н‚[а,Б]) = со сВеБ(0,2(х)) +...+ с„ сВеБ(и,2(Х)), со,..., с, Е Ающт. 
Метод 2. 


1. По методу Галеркина в пространстве [2(а,Ь;р) преобразовать ЛИУМК Ву + 
+ Ну]| + /= 0 системы (6) в многочлен 


А, = Саейет_ Го(а, Б;р)(В у + НУ] + = 0, [а,6], п). (22) 


2. Преобразовать уравнение А у + [у] + = 0 системы (6) в уравнение вида (8) с 
многочленом А„ (22) 
Ану + Цу|+8=0. (23) 


3. По методу Галеркина в пространстве /2(а,Б;р) преобразовать уравнение (23) 
в искомый многочлен 


у, = дает _ [2(а,6;р) (Аку + Цу|+2=0, [а,6], п). (24) 


6 Алгоритм — реализация метода 2 


ПМРОТ = (Ау + Цу] + & =0, А = зо уе(В у + Му] + = 0), [а,6], п). 
Алгоритм 2. 
1. Преобразовать ЛИУМК Ву + Н\у| += 0 системы (6) в многочлен 


Аи = зоуе(5 В (у„ е Ниа,Б]) + Н [(”, е Нща,6])] +1=0). (25) 
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2. Преобразовать уравнение А у + Цу] + © = 0 системы (6) в уравнение вида (8) с 
многочленом А)» (25) 
Ау + Цу| +#=0. (26) 


3. Преобразовать ЛИУМК (26) в многочлен 


а|еотийт_2(А у + Цу]| + 2 =0, А = зоуе(Ву + Му] += 0), [а,6], п) = 
ООТРОТ = у, = $91%е(5,„[ А» (у, Е Ни[а,6]) + [ (», Е На,6])] + 2] = 0). (27) 


Лемма 3. Метод 2 эквивалентен алгоритму 2 — тождественны многочлены 
(21) и (27). 

Доказательство. Согласно определению метода Галеркина в пространстве /2 (а,Б;р), 
многочлен А» (25) тождественен многочлену А) (22). Поэтому тождественны уравнения 
(23) и (26). Следовательно, согласно определению метода Галеркина в пространстве 
[2(а,6;р), многочлен у, (21) тождественен многочлену у/„ (27). 

Оптимальность алгоритма 2. В пространстве Г2(а,6;р) оператор 5» (21) — 
наилучший аппарат для аппроксимации функций (у(х), х Ее [а,Б]) е 1[Г.(а,6:р) 
алгебраическими многочленами — справедливы тождества 


[у — $=[У] [22,5;р) = 1 о,.. св || У(%) — (Со... сих”) [жару |5 ежа, = 1. 


Определение 1. Система уравнений (6) принадлежит классу А'5 . <=> 
Коэффициент оптимальности (3) варианта (21) метода Галеркина в простран- 
стве [2(а,Б;р) на этой системе уравнений в пространстве Х = Си] ограничен 


С„(Саегкт, А у + Цу| + ® =0, А = зоуе(В у + Ну| += 0), Х) | уб® - у ||х/ 
ШЕ 0... са || У - (со +... сих”) [к= 0@). 


В пространстве Сри,ь оператор 5„ — оптимальный среди операторов проектиро- 
вания функций (у(х), х Е [а,6]) е Сщы в пространство алгебраических многочленов и 
справедливо тождество 


|5» [саы] = (4/<^) ни) + ОСТ), (1) < 3. 
Поэтому класс А'; систем уравнений вида (6) достаточно широкий. 
Класс А'5 оценивают условия основной теоремы сходимости метода Галеркина (21). 
7 Оптимальность метода 1 
Лемма 4. Аппроксимация ЛИУМКАу-+ Цу]+8=0 (4) по методу [3] 


А (,ЕР,) + Цул Е Р,] ++ (Ешл Е Нил [а,Ь]) = 0, (28) 
где т = 4е2(А (у, ЕР.) +ЦШу, Е Р,„|+8) (29) 
(»ЕР,) = с +с1х+...+ сих", со, с1,,..., © Е Аютщ, (30) 


(Етл Е Нти [а,6]) = <! сВеБ (и-+1,2(х)) -...+ ти-„ свеб(т,2(х)), т! ,..., ти Е Аш, (31) 
эквивалентна системе: 
— аппроксимация этого ЛИУМК по методу Галеркина в пространстве [2(а,Б;р) 


5[А (у, е Наа,Б]) + Ц (» е Нша,6])] + 8] = 0, (32) 

— система уравнений 
{акА Уп + 7Й] Уп ] + 2 [а,6] ) + Тм = 0 } Е (33) 
где у = зО1уе(5„[А (у, Е Ни[а,6]) + ЦО» Е Ниа,6])] + ] = 0) (34) 
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Доказательство. Для ЛИУМК (4) левая часть уравнения (28) — многочлен порядка 
т (29). Многочлен Р„, — порядка т тождественен нулю тогда и только тогда, когда нулю 
тождественны его коэффициенты Фурье — Чебышева порядка 


1=0,...,т на отрезке [а,6] — а{ Ри, [а,6]) = 0, 1=0,...‚т. 


Следовательно, уравнение (28) эквивалентно системе уравнений 


{аКА (у е Р»„) + МЦуле Р,] +5 + (Етл Е Ни [а,6]), [а,6] ) = 0}:=0" (35) 
относительно коэффициентов многочленов 
(и ЕР,,) (30) и (Ешл Е Ним [а,6]) (36) 


Функционал а{Ри, [а,6]) — линейный и для него справедливо тождество 
а(со свеб(0,2(х)) +...+ с; свеб(т,2(х)) +...+ с„ свеБ(и,2(%)), [а,6]) = с;. 


Следовательно, справедливы тождества 
ак(Етл Е Ним [а,6]), [а,6]) Е 0, 1= 0: п, 
ак((Етл Е Нти [а,6]), [а,6]) = лт-ь = И+1,..., Т. 
Поэтому первые и + 1 уравнения системы (35) имеют вид 


а(А Е Р,) +Цу ЕР, | + в, [а,6]) = 0,1=0,...., п. (37) 


Элементы пространства Н„[а,6] — алгебраические многочлены порядка п и толь- 
ко эти многочлены. Поэтому многочлен у, (34) тождественен многочлену 


уи = $01Уе( 5„[ А (ЕР) +ЦуЕР, | +8] =0), 5, : Е2(а,б; р) —> Ниа,6]. (38) 


Следовательно: 

— уравнение (37) эквивалентно системе уравнений (36), 

— решение системы (36) — коэффициенты многочлена у» (37), 

— результат подстановки многочлена у, (37) в систему уравнений (35) вместо 
У, ЕР») — система уравнений (33) ип + 1 тождество 0 = 0. 

Лемма 5. //усть уравнение Ау + Цу]-+8=0 (4) - ЛИУМК. 

Тогда на этом уравнении алгоритм [3] эквивалентен методу Галеркина в 
пространстве Г2(а,Б;р) — тождественны многочлены у, (34) и 

у, = тефо4з] (А у + Ду] +&=0, [а,6], п). (39) 

Доказательство. Согласно лемме 4, уравнение (28) эквивалентно системе: урав- 
нение (32) и система уравнений (33). 

Уравнение (32) является аппроксимацией уравнения (4) по методу Галеркина в 
пространстве Г2(а,5;р). 

Уравнение / = 1,...,-и системы (33) определяет коэффициент т; через коэффи- 
циенты многочлена у, (34). Следовательно: 

— уравнения (33) являются тождествами, 

— многочлен у, (38) тождествен решению у, (34) уравнения (32). 

Лемма 6. Алгоритм 1 эквивалентен алгоритму 2 — тождественны многочлены 
У, (12) и (27). 

Доказательство. Согласно лемме 5, тождественны многочлены А) (25) и (10). 

Поэтому тождественны уравнения (11) и (26). 

Следовательно, согласно лемме 5, тождественны многочлены у» (12) и (27). 

Теорема 5. Метод 1 эквивалентен методу 2 — тождественны многочлены уп 
(9) и (21). 

Доказательство. Согласно лемме 1, тождественны многочлены у, (9) и (12). 
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Согласно лемме 6, тождественны многочлены у/ (12) и (27). 

Согласно лемме 3, тождественны многочлены у/ (27) и (21). 

Теорема 6. Пусть система уравнений (6) принадлежит классу А'5 . 

Тогда коэффициент оптимальности (3) алгоритма [ на этой системе уравне- 
ний в пространстве Х = Си, ограничен. 


С; аюогИт_1, Ау + Ду] + = 0, А = зоуе( Ву+ Ну] +Х=0), Х) = 
|У— у» |х/ Ш о, ...с п || У(®) — (со +... си х')[х=О (1). 


Доказательство. Согласно лемме 6, тождественны многочлены у, (12) и (27). 

Согласно лемме 3, тождественны многочлены у, (27) и (21). По определению 1 
класса А'; ‚ коэффициент оптимальности метода Галеркина в пространстве /2(а,5;р) 
(21) на системе уравнений этого класса в пространстве Су] ограничен. 


Заключение 


Построенная АРГАМ-процедура доказывает эффективность представленных 
инструментальных средств для создания процедуры СКА для решения интегральных 
уравнений отдельного типа. 

Эта АРГАМ-процедура в случае системы интегральных уравнений вида (6) ре- 
шает задачу Дзядыка [2]. На основании а-метода построить эффективные алгоритмы 
для вычисления алгебраических многочленов — аппроксимации решения у(х), х е [а.Ъ] 
функциональных уравнений оптимальной в пространстве Саь|. 


Дополнение 


1. Алгоритм 1". Модификация алгоритма 1. В алгоритме [3] базис пространства 
[2(а,Б;р) заменен на базис пространства Гильберта Наы. 

2. Для алгоритма 1" имеют место аналоги теорем 5 и 6. В этих аналогах прост- 
ранство [2(а,Ъ; р) заменено на пространство та,ы. 

3. Модификация АРГАМ-процедуры. Реализация алгоритма 1". Реализация ал- 
горитма [3] заменена реализацией модификации этого алгоритма с базисом простран- 
ства Н! [а,5] . 

4. Распространение алгоритма 1 на интегральные уравнения вида (6), где коэф- 
фициенты — решения ЛИУМК (4) или многочлены, эквивалентно соответствующему 
варианту метода Галеркина в пространстве Г2(а,Б;р). 

5. Этими инструментальными средствами мы построили процедуры для ре- 
шения интегральных уравнений вида 


Ау + Цу] + 8 =ХДх, Шу) иАу + Цу| + &=НК_, У}, 


где Г, Н — линейные интегральные операторы вида (5) с многочленными коэффи- 
циентами и нелинейность — Хх, у) — функция или многочлен. 

Эти процедуры вычисляют многочлен у„ — аппроксимацию решения исходного 
уравнения оптимальную (3) в пространстве Са. 
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П.М. Денисенко 

Розвязування 1нтегральних р!внянь в системах комп’ютерно! алгебри 

Представлено 1нструментальн! засоби для системи алгебрайчного програмування АР$, Маре, Мафетайса 
та 1нших систем комп’ютерно] алгебри. Цими засобами конструюють комп’ютерн! програми для таких 
систем. Ц! програми мають наступн1 вх та вихд: 


МРОТ = (И у]=0, [4,6], п), где Ну] =Хх, бо), [| Кб, у@)) а), 
ОПТРОТ = у, = с с,х+ ...+ сх" = у(х ) = зо№е( Е[у]=0), хе [а,Б] , 
|| УС%) - у.) [атаьы / ПаВе о, еп || УСО - (С%+... + сих”) [свы = О. 


Цими засобами сконструйована процедура системи АР$. Ця процедура розв’язуе лиийн! 1нтегральн! 
равняння виду 


Ау+Цу]+8=0, де А = зоуе(В у - НУ] + /=0), Ду] =Киу] +...+ КЫУ], 
НОУ] = Кн +... + КБ, КО = [у КК У, 
ядра КХх,1) и коефипенти сх), Хх), ©, В, /— багаточлени. Доведена ефективнасть процедури. Доведена 


ефективнисть алгоритму та методу — основ процедури. Ця процедура доводить ефективнисть побудова- 
них засоблв. 


Р.М№. Рет5епКо 

Тре Пцерта! Едиайоп$ Зоушо ш бе Сотршег А1хеБга Зу$етб$ 

\'е ргезещед фе еЁйслепе таетайса!| аррагайаз (Фе 10015) Рог Фе А]ебга!с Ргоэтатите Зузет (АРЗ), 
Маре, Мафетайса ап4 оег сотрщег а1себга зузетз. ТБ15 {0015 сопзгасе Бе сотршег ргозгаплз. ТВезе 
ргоетатаз Бауе: 


МРОТ = (И у] = 0, [4,6] п), где Ну] =ЛХх, уб%), |? Кб, У@)) 4), 
ОПТРИТ = у, = с +с,х+ ...+ с, х" = у(х ) = з9уе(Ё [У] = 0),х Е [а,Ь], 
|| 5%) - у,б® Даы / На, о, це и [У - (соч... + с, | свы= О). 


Озше Фезе 10015 ме сопзбаце4 Фе ргоседиге Юг АР$. 'ТЬ15 ргоседиге 5о[уез Фе Ппеаг ицеога| едиайоп$ 


Ау+Цу]+8=0, \Пеге 4 = зо№е(В у + Ну ] += 0), Шу] = КЦу] +...+ КУ], 
НУ] = Кн + ... КБ, КО = [у КК У, 
Фе Кегп$ К (х,0) ап4 Фе соеЁйслеп сх), 4х), ©, В ‚ Г- Фе роПпопиа15. \е ргоуе4 е еЁЙсепсу оЁ 15 


ргоседиге. \!е ргоуе4 пе еЁйслепсу оЁ 1$ те#фоа Ус 1$ Фе ЪБа$15 ог 15 ргоседите. 'ТЬ$ ргоседиге ргоуез 
Фе еЁЙстепсу оЁ Шезе 10015. 


Статья поступила в редакцию 14.06.2010. 
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